NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 



TECHNICAL NOTE 2431 


SKIN FRICTION OF INCOMPRESSIBLE TURBULENT BOUNDARY 

LAYERS UNDER ADVERSE PRESSURE GRADIENTS 

By Fabio R. Goldschmied 

Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 



ATWOOD 


2192 


TECH LIBRARY KAFB, NBA 



DbSbta 


1 NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2431 


SKIN FRICTION OF INCOMPRESSIBLE TURBULENT BOUNDARY 
LAYERS UNDER ADVERSE PRESSURE GRADIENTS 
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SUMMARY 

Experimental data for skin friction of turbulent boundary layers 
under adverse pressure gradients from several sources are presented in 
graphical form. Data obtained by the momentum-balance method are shown 
to follow a trend opposite to that of data obtained by hot-wire and 
heat-transfer methods. The momentum equation and total-pressure- 
measuring techniques are discussed in relation to skin-friction 
calciilations . 

The conclusion is reached that momentum method is of questionable 
value in regions of adverse gradient and should not be relied upon 
because the sensitivity to small measuring errors and to deviations from 
assumed flow conditions in test channels is too great. A new integral 
energy parameter is introduced and its relation to skin-friction data 
is demonstrated on the basis of available material. A new momentum 
thickness which includes the turbulent momentum contribution, is also 
introduced, and it is shown that this thickness may be as much as 

T~ percent lower than the conventional momentum thickness, near the 
turbulent separation point. 


INTRODUCTION 

The present knowledge concerning the development of turbulent 
boundary layers under adverse pressure gradients and approaching 
separation is far from adequate. In particular, the skin friction 
under adverse pressure gradients cannot be determined with reasonable 
accuracy because of the lack of enough reliable experimental and 
theoretical information. Only through positive boundary- layer- control 
devices, such as suction slots, suction areas, vortex generators, and 
so forth, can the design engineer be reasonably assured that the boundary 
layer will develop in a desired manner. 

The present analysis was prepared at the NACA Lewis laboratory to 
assemble and to review in a critical manner the available information 
concerning the skin- friction trend of an incompressible turbulent 
boundary layer in regions of adverse pressure gradients in order to 
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show the conflicting data obtained, by different experimental methods 
and to bring out the important characteristics and parameters affecting 
skin- friction development . 

Because of the absence of a known fundamental relation between the 
mean velocity gradient and the statistical quantities describing the 
turbulent field in the general case of shear flow,, the theoretical 
approach to the study of turbulent boundary layers has been based almost 
exclusively on integral methodsj possible solutions have been sought that 
would satisfy, not the local differential relations as given by the 
Reynolds equations (reference l) , but only "over-all averages" of momen- 
tum (reference 2) and energy (references 3 and 4) across a boundary- layer 
section. 

In general, the averaging across the boundary- layer section may be 
weighted by multiplying the equation of motion before integration by 
some arbitrary function chosen to give greatest importance to some 
particular area of the bound ary- layer section, such as the area closest 
to the wall or the area nearest to the free stream. Such weighted 
integral equations are developed in references 5 to 7. 

These methods provide one or more algebraic relations between wall 
friction and an equal number of parameters employed to describe the mean 
velocity profile. Thus, even for cases of zero pressure gradient, a 
solution for boundary- layer development requires an empirical relation 
for skin friction. 

A sufficient quantity of data available for skin friction of low- 
speed turbulent boundary layers under zero pressure gradients (flat- 
plate case) has enabled several authors to develop empirical formulas 
for skin friction as a function of the Reynolds number based on momentum 
thickness. These formulas are strictly applicable only to the flat-plate 
case. German authors generally use the Schultz-Grunow formula (refer- 
ence 8 ) j British and American authors, the- Squire and Young (reference 9) 
and the Falkner (reference 10 ) formulas. 

When the turbulent boundary layer develops under an adverse pressure 
gradient, both skin friction and mean velocity profile depend in some 
manner on pressure gradient. In laminar flow, the pressure gradient 
enters the analysis through a single d im ensionless group called the 
Pohlhausen \ (reference ll) . Doenhoff and Tetervin (reference 12), Burl 
(reference 13), Garner (reference 14), and Kalikhman (reference 15) have 
proposed the use of various quantities, all analogous to the Po hlhau sen X, 
to describe the effect of pressure gradient on the turbulent boundary 
layer. 


The following experimental techniques have been enployed for the 
measurement of turbulent skin friction in an adverse pressure gradient: 
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(a) Total-head survey: Mean velocity profiles, measured hy total- 

head surveys, permit solving the integrated momentum equation for skin 
friction. This method is in principle simply the reverse of the 
theoretical approach discussed previously, -whereby a skin-friction law 
is assumed in order to compute mean profiles and development. 

(b) Heat transfer (reference 16) : A small heating element is 

inlaid flush with the test surface. Heat is convected away from this 
element by the flow in the laminar sublayer next to the wall. The film 
coefficient may be measured and related by calibration to skin friction. 

(c) Hot-wire anemometer (references 17 and 18): Velocity fluctu- 

ations affect, by convective cooling, the electric resistivity of a small 
heated wire. Turbulent fluctuations in the boundary layer may therefore 
be measured to provide an evaluation of the Reynolds stress close to the 
test surface. 

This report is primarily concerned with a comparative discussion of 
these techniques, the data obtained by the various investigators, and the 
resulting theoretical implications. 


DISCUSSION 

Influence of Surface Roughness, Free-stream Turbulence 
and Transition History 

It is well known that skin-friction values depend on surface- 
roughness conditions (reference 19 to 21); consideration must also be 
given, however, to the influence of free-stream turbulence on skin- 
friction development for a given chordwise position of the transition 
from laminar to turbulent flow. The transition history of a turbulent 
boundary layer has been investigated experimentally (reference 22) and 
it has been found, in the case of zero pressure gradient, that for a 
given Reynolds number R^, the state of the boundary layer is the same, 
with respect to mean velocity profile and energy spectrum shape, regard- 
less of the m a nn er in which transition has been caused, provided that a 
sufficient length of run is allowed for settling. This investigation of 
transition-history effects has proved to be very important; because here- 
tofore it had not been apparent that the turbulent boundary- layer develop- 
ment is independent of the type of transition, which may occur in a 
multitude of different ways. 

Surface roughness and free-stream turbulence both influence the 
energy transfer from mean to turbulent flow by increasing the turbulent 
shear and modifying the mean velocity gradient. An example of the 
influence of free-stream turbulence on the skin friction of a turbulent 
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"boundary layer un der adverse pressure gradient, as presented "by Tillman 
in reference 21, is shown in figure 1. Percentage wise, the difference 
in skin-friction values between runs with and without a turbulence grid 
is quite appreciable and comparison between results of different authors 
cannot be properly made unless turbulence conditions are exactly speci- 
fied by relevant statistical quantities along the flow path, which has 
very seldom been done in the past. 

The difference in mean velocity profiles at the same station for a 
boundary layer on a flat plate with and without a turbulence grid, as 
presented by Wieghardt (reference 20) is shown in figure 2» The boundary 
layer thickness is respectively 80 and 50 millimeters, as given in ref- 
erence 20. Consequently de/dx may also markedly depend on free- stream 
turbulence. (Symbols are defined in appendix A.) 

Most data have been obtained under low- turbulence conditions and 
great caution must be used in applications to high-turbulence conditions, 
such as exist in turbo -machinery. 


Skin Friction in Adverse Pressure Gradient by Momentum Method 

A number of authors, among them Gruschwitz (reference 23) , Tillman 
"(reference 21), Wieghardt (reference 24), and Wieghardt and Tillman 
(reference 25) have measured the skin friction of a turbulent boundary 
layer under adverse pressure gradient by the Karman momentum method. 

The Karmen momentum equation may be written: 


s + e v (E+2) - °f 0 (1) 

As shown in appendix B, this equation may be obtained by integrating the 
Reynolds differential equation through the boundary layer, subject to the 
assumptions arising from Prandtl's concept of a thin boundary layer as 
applied to turbulent motion. Thus, the mean pressure is taken to be a 
function of x alone, and the x-derivatives of both mean or fluctuating 
quantities are neglected in comparison with the y-derivatives of the 
same quantities. This second assumption is doubtless valid for mean flow 
q uant ities, but for fluctuating quantities justification must be based on 
analysis of experimental evidence. An analysis of this sort will be pre- 
sented subsequently. 

Total -pressure surveys and wall static-pressure measurements suffice 
to determine the left-hand side of equation (l), which is then identified 

with the wall skin- friction coefficient C-p . 

r 0 
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For illustration, curves of the skin-friction parameter 



against x are plotted in figures 3(a) and 3(h) from data of refer- 
ences 25 and 18, respectively. Another example of the variation of 
Cfo with x is given in figure 1 from data of reference 21. For the 
cases of figure 3 which represent different pressure gradients, neither 
surface conditions nor free-stream turbulence have been specified, but 
transition has been artificially forced, always at the same position, by 
means of a wire set in place immediately after a suction slot. Pressure 
distributions for all cases illustrated can be found in the references 
noted. Faliner ' s curve has been shown in figure 3 for convenient com- 
parison with what would be corresponding flat-plate values. The curves 
all show a considerable increase in the direction of flow, whereas by 
laminar analogy it was expected that the skin friction would decrease 
monotonically to zero at separation. 


The apparent increase of skin friction in the flow direction, 
which at first appears to be an established experimental fact, has 
fascinated many researchers. Doenhoff and Tetervin (reference 12) 
recognized the contrast and the possibility that the total-pressure 
measurements were too high because of relatively high turbulent fluctu- 
ations, but nevertheless accepted the increasing trend as substantially 
correct. Garner (reference 14), who rejected the skin- friction data of 
reference 16, also accepted the increasing skin-friction trend. 

Gruschwitz (reference 23) and Tillman (reference 21) looked for possible 
causes of error. 

No practical use, however, has been made of these data for appli- 
cation to solutions of the K^rm^n momentum equation for 9; only flat- 
plate empirical formulas based on Rq have been used for this purpose, 

in cases of pressure gradient, on the assumption that 9 does not 
depend strongly on the precise form of the function Cf^(x). Doenhoff 

and Tetervin (reference 12) and Kalikhman (reference 15) use the Squire 
and Young formula (reference 9). Garner (reference 14) uses the Falkner 
formula (reference 10) while the Schultz-Grunow (reference 8) formula 
is used by German authors. 

Verification that the momentum thickness 9 computed from the 
K^rm^n equation is insensitive to C-, may be obtained from the data 

r 0 

of Schubauer (reference 18 )j the momentum thickness 9 has been computed 
with his experimental values of H, by using several sets of values for 
Cf of which two are of interest here: 

(1) The experimental values of Cfg computed according to the 
K^rm^n momentum equation from the faired experimental 9 points 

(2) Falkner ' s skin-friction values based only on Rq and with 
pressure gradient effects neglected 
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These two values of (together -with a curve obtained with a 

hot-wire anemometer to be discussed later) are plotted in figure 4. 

The difference in the resulting 9 values, shown in figure 5, is 
proportionately quite small, so that it seems reasonable to use flat- 
plate skin-friction formulas, such as Falkner ' s when 9 is to be 
computed. In other words, at least in the range of these data, 9 
is really insensitive to changes of C-^ . (in fig. 5 a set of points 

of "turbulent momentum t hickn ess" and a set of 9 values computed on the 
basis of hot-wire skin friction are also presented for later discussion.) 
On the other hand, when skin- friction values are computed from experi- 
mental 9 data by using the K^rmdin momentum equation as an algebraic 
equation for it appears that Cf Q is very sensitive to small 

changes in do/cbc. 


is a universal function of 


One-parameter correlation for C-p ^. - From laminar analogy, several 
authors have attempted to show a one-parameter correlation for Cf^ as 
obtained by the momentum method. Gamer (reference 14) attempted to 
show that ~ C-p 

6 f 0 « . 
inferred from some experimental data that — Eg ' is constant and 

equal to 0.007623 for ORq 1 / 6 * -0.01. The data of references 25 and 

18, plotted in the manner of Garner, are shown in figure 6 and seem 
clearly to disprove his contention. 


GRq 1 / 6 21 and 


Wieghardt and Tillman (reference 25) and Tillman (reference 21) 

U' 

attempted to correlate Cf^ with the pressure parameter 9 — . Flg- 

U' 

ure 7 shows curves of C-p against 9 — , from data of reference 21. 


L 0 


U 


1/6 


against a 


Figure 8 shows curves of skin-friction parameter Cj=^ 

similar pressure parameter ijr 21 from the data of references 25 and 

U 

18 ) it appears that a typical curve Is one with a return loop, so that, 
for a given value of the pressure parameter, two values of the skin- 
friction parameter are possible. Furthermore, the points are scattered 
over a wide region instead of forming a universal curve. 


These data indicate most clearly that the skin- friction parameter 
does not depend on the local values of the pressure parameter in the 
manner that the skin friction depends on local values of Rq in the case 
of zero pressure gradient. 
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Two-parameter correlation of C ^. - For a two-parameter correlation 

of Cfgj a new second parameter must be selected as none appears in the 

literature. Surface roughness and free-stream turbulence have an effect 
on the boundary- layer development because the rate of energy transfer 
from mean to turbulent flow is modified across the boundary- layer section. 
It was therefore believed reasonable to select a parameter indicative of 
the profile of energy transfer, in dimensionless form, such as the right- 
hand side of the mean energy equation (derived in appendix C). This para- 
meter, denoted by E, is 



where the energy deficiency thickness i)r t 

P8 


% 


u 

U 


0 


1 - 


u^ 
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As was expected, the skin-friction parameter CL. R did not 

J 0 ® 

correlate against local values of the parameter E. If, however. Of 


0 


is assumed to be a function of R , — i|r, and E, then a plot of CL. R^ 1 / 6 

y U r o ® 


against E for constant values of 
degree of correlation. 


IT 

U 


■>|r might be expected to show some 


Typical results of such an attempt at two-parameter correlation are 
shown in figure 9. It would appear that no matter how small the pres- 
sure parameter is chosen RqV® increases monotonically as a function 

of E, possibly along the straight lines shown in figure 9. An appreci- 
able amount of scatter of the points occurs about these lines, but the 
scatter is not too large to be attributed to experimental or computational 
error. 
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The foregoing approach is of value only insofar as the momentum 
measurements hy total-head surveys are substantially correct and the 
quantity heretofore denoted Cj^ maintains a physical significance as 

the right-hand side of the momentum equation. 


Hot-Wire and Heat-Transfer Measurements 


Experimental methods have "been developed that allow evaluation Of 
the turbulent shear stress in the vicinity of the wall in regions of 
adverse pressure gradient. This turbulent shear is considered to be 
transmitted across the laminar sublayer to the wall according to the 

relation =* obtained from Reynold’s equation and evaluated at 

the wall. Schubauer and Klebanoff (reference 18) have shown that the 
skin friction, as determined by an X-type hot-wire set as close as 
possible to the wall, according to the scheme just described, exhibits 
a monotonic decline in the region of adverse pressure gradient down to 
zero at separation. 


Ludwieg (reference 16) recently developed a heat-transfer type 
instrument, which, when calibrated dynamically, is capable of direct 
measurement of skin-friction magnitude and direction. A small heating 
element, flush with the wall is cooled by convection, losing heat 
through the laminar sublayer. The heat-transfer coefficient is 
evaluated experimentally and related by calibration to the local skin- 
friction coefficient. 


This method of measurement, though requiring careful calibration, 
would appear capable of providing a direct measure of the wall shear- 
ing stress and has the additional important advantage that its use 
probably involves negligible aerodynamic disturbance of the flow, since 
the instrument is flush with the wall and the velocity profile is 
independent of the temperature profile for substantially constant prop- 
erty values of the fluid. 


Certain data have been obtained with this instrument by Ludwieg 
and Tillman (reference 26) that show a monotonic ally declining skin- 
friction trend for two cases of adverse pressure gradients. These 
data cannot be regarded as conclusive since pressure gradients were not 
specified except as to sign, boundary- layer thicknesses were not given, 
and the separation points were not indicated. 

In figure 10, the data from reference 18 and reference 26 are 
plotted against R^j for comparison, the empirical flat-plate skin- 

friction formulas of references 8. to 10 are also represented. The hot- 
wire points lie on a curve which seems to be too high to assure a smooth 
joining with the flat-plate curves. On the other hand, the data of ref- 
erence 26 join smoothly with the flat-plate curves. 
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PROBLEM OF CONFLICTING SKIN-FRICTION EXULTS 

Inasmuch as the momentum, hot-wire, and heat-transfer data dis- 
cussed in the preceding sections are presented "by different authors, 
such data are not necessarily coup arable or applicable to the same 
phenomenon. In Schubauer's case, however, it is possible to compute 
the momentum skin friction from the data available and to compare it 
with the hot-wire results} the data points are plotted in figure 4, and 
show that a rising momentum-skin-friction trend is concomitant with a 
decreasing hot-wire skin-friction trend. Such a check could not be 
made from the Ludwieg-Tillman data because of insufficient information. 

This discussion suffices to show that the experimental evidence 
concerning the skin friction in an adverse gradient is glaringly con- 
tradictory. The trend shown by the hot-wire method (and supported by 
the heat-transfer data) may be correct} the trend shown by the momentum 
survey method may be correct; or both may be in error. The following 
paragraphs present a preliminary study of the errors inherent in the 
current methods of obtaining skin-friction results in an adverse pres- 
sure gradient. 

Possible inadequacies of momentum theory based on two-dimensional 
stationary mean flow. - The K&rm&n momentum equation may not be 
adequate for adverse-pressure-gradient conditions. What has been 
computed as skin- friction coefficient may actually be the sum of 
several integral terms, or in other words, the well-known Prandtl 
boundary- layer assumptions may not satisfactorily apply. No experi- 
mental or theoretical evidence has been presented that the Karmen 
momentum equation does provide an adequate momentum balance for the 
turbulent boundary layer under adverse pressure' gradient and approaching 
separation, whereas the hot-wire and heat-transfer data tend to show 
that it does not. 

It is therefore necessary to return to the complete equations of 
motion for two-dimensional mean flow with three-dimensional turbulence, 
with the assumption that mean turbulent quantities do not vary in the 
z-direction. In appendix B, the following integral momentum condition 
is derived from the equations previously described; 


^t U 1 - 

dx“ + IT 


9 t (Ht+2) 


C-p + A + B 
r 0 


( 2 ) 


e. 
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U 
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The conventional laminar momentum- thickness expression is modified 
to include turhulent-momentum contributions,, -which also means that the 
shape factor is modified accordingly because 5* remains unchanged. 

(The turbulent fluctuations by definition do not contribute to mass flov.) 


In isotropic turbulence, -where a T yy> and in zero turbulence, 

where 5 Tyy = 0, 0 » 9^.; however, it is well known that the turbulent 

field must be strongly nonisotropic in shear flov. Only at the free- 


stream boundary and at the edge of the laminar sublayer will 
actually be zero. 


xx 


yy 


In figures 11 and 12 the mean and turbulent profiles, from refer- 
ences 18 and 20 are shown. Values of 0 and 0^ computed from these 

data are presented in the following table for. data from reference 20: 


Turbulence grid 

0 

e t 

0 - 0 t 

y i nn 

R 

X 

Fig. 

St 

(percent) 



With 

4.48 

4.30 

4.18 

6.25 X 10 6 

l?(a) 

Without 

5.39 

5.22 

3.26 

6.25 X 10 6 

12(b) 
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and values "based on reference 18 are presented herein: 


X 

9 

e t 

e - e + 

—5 X 100 

w t 

(percent) 

s e 

Fig. 

17.5 

0.2318 

0.2202 

5.27 

0.213 X 10 5 

11(a) 

20.0 

.3302 

.3146 

4.96 

.281 

11(b) 

21.0 

.3984 

.3796 

4.96 

.321 

11(c) 

22.5 

.6308 

.5988 

5.34 

.472 

11(d) 

23.5 

.7228 

.6916 

4.52 

.511 

11(e) 

24.5 

.9468 

.9076 

4.32 

.634 

11(f) 

25.4 

1.1812 

1.098 

7.57 

.766 

11(g) 


From these tabulations it is apparent that the turbulent momentum 
thickness 9^ can be over 7.5 percent smaller than 9) the difference 

between dB^/cbc and d9/dx may be larger , but it is impossible to 


estimate how much. In figure 5, 9-^ is compared with the experimental 
values of 9 and with values of 9 conputed from the formula 

Cf Rq 1 / 6 a 0.006534 (Falkner's formula), and from hot-wire results. 


On the right side of the conplete momentum equation there are then 
three terms, of which the first is the skin- friction coefficient C^, 

the second is the contribution from deviation of the turbulent flow- 
field from the Prandtl assumptions A, and the third is the contribution 
from the deviation of the mean flow field from the Prandtl assumptions B. 
A measurable pressure gradient must exist across the boundary- layer 
section if terms A, B, or the sum assume an appreciable magnitude. In 
other words, the transverse pressure gradient can originate from 
longitudinal gradients of the turbulent flow field, the mean flow field, 
or both. 


According to the Prandtl approach, longitudinal gradients can be 
neglected in comparison with transverse gradients. Prandtl, however, 
referred to laminar flow fields; in turbulent regimes, his assumptions 
can be expected to apply only to the mean flow field, as Liepmann and 
Laufer (reference 26) point out. 

On the basis of direct pressure measurements, Schubauer and 
KLebanoff (reference 18) state that only "barely detectable" pressure 
differences were found across the boundary layer at some stations; it 
therefore seems safe to conclude that no appreciable transverse gradient 
has been found experimentally. On the other hand, computations of terms 
A and B should result in agreement with this conclusion; that is, A + B 
should be negligible in comparison with Cp . 
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On the "basis of the extensive experimental material presented in 
reference 18, this computation of terms A and B has been attempted. 
Because of the difficulties of obtaining first and second x-derivatives 
graphically, however, the accuracy of the results is not sufficient to 
allow their presentation; A and B are probably indeed negligible up 
to the vicinity of the separation point. If A and B are negligible, 
the right side of the momentum equation is represented by Cf^ alone, 
at least for engineering purposes, and the discrepancies between momen- 
tum skin-friction data and hot-wire and heat-transfer skin-friction can- 
not be attributed to inadequate momentum theory for two-dimensional 
flow. 


As was mentioned in the preceeding section, hot-wire measurement 
of wall friction involves, according to theoretical considerations, a 
straight-line extrapolation from a point on the experimentally determined 
Reynolds stress profile. This point is taken as the point of tangency 

^ T xy dp 

of the Reynolds stress profile to the straight line of slope — » —• 

oy dx 

Thus it is assumed that down to the immediate vicinity of the wall, the 
Reynolds shear is much larger than the molecular shear, and that in this 
•- vicinity the curvature of the total shear profile is negligible. The 
accuracy of these assumptions remains to be evaluated experimentally. 


Inadequacy of experimental technique : Velocity measurements in a 

turbulent field can be made either by a total-head tube or by a hot- 
wire anemometer. Liepmann and Laufer (reference 27) show that the 
total-pressure-probe value is consistently higher than the hot-wire 
value by an appreciable amount. If a total -pressure tube is con- 
sidered as reading the sum of the mean velocity and of the root- 
mean-square fluctuating velocity, the effect of this error on the 
determination of 8 will depend, among other things, on the 
relative fluctuation profile; that is, on the profile of u'/u. It is 
well known that u'/n can assume values up to 0.40 in the vicinity of 
the wall. Inasmuch as d9/dx Is wanted rather than 0 itself, it is 
not clear what contribution to a possible C r error Is given by high 

x 0 

total-pressure readings. 


Possibility of secondary flows or violent unsteadiness : As the 

separation region is approached, secondary flows may be induced in the 
boundary layer, which invalidates the two-dimensional assumptions. 
Tillman made an investigation of this effect in rectangular test channels 
for his doctoral thesis (Gottingen, 1947) and states that an appraisal of 
the effect of secondary flows gives a % value lower by some 40 per- 
cent, as discussed in reference 26. There is a possibility that the 
turbulent fluctuations become large near separation; such large fluctu- 
ations might impart a violently unsteady character to the flow and cause 
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large cross flows. The decidedly nonlinear "behavior of turbulent flow 
is evidenced by the recent observations (reference 22), that turbulent 
fluctuations occur actually in gusts or bursts in turbulent shear flows 
with a free boundary. In this event, none of the methods discussed 
would be adequate in the vicinity of separation. 

This discussion of the methods for determining skin friction of 
turbulent boundary layer under adverse pressure gradients and approach- 
ing separation thus indicates that the momentum method is of question- 
able value because it is too sensitive to errors in the determination 
of 9 values and because it assumes two- dimensionality of flow, which 
may be difficult to achieve in test channels. It would therefore seem 
plausible, on the basis of the evidence available, to consider the skin- 
friction trends observed by the hot-wire and heat-transfer methods to 
be more reliable than the momentum trend. The evidence should not, 
however, be regarded as conclusive. 


RELATION OF SKIN FRICTION TO RATIO OF ENERGY QUANTITIES AND TO 
MEAN VELOCITY-PROFILE PARAMETER 


For the data of reference 18, a certain integral-energy ratio shows 
exactly the same trend as the hot-wire skin friction, which may be coin- 
cidental inasmuch as no supporting evidence as yet exists. 


The quantity 


r 8 


i pu 3 dy 


nx p6 


JO JO 


Vs! ■v 


1 


is plotted against the experimental hot-wire 



in figure 13. 


The similarity of the curves is close over the entire range, so 
that the points bracket the 45° line quite well; furthermore, this 
similarity also extends into the favorable-gradient region (not shown 
in fig. 13). In terms of boundary-layer thicknesses, this similarity 

indicates that C„ is proportional to - Jl- - 2~1 . Thus, in 

f o [t * . J 

physical terms, the ratio of the volume integral of the turbulent energy 
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production to the section integral of the kinetic energy appears to he 
related to the local skin friction as determined hy the hot-wire method. 

From their heat-transfer data, Ludwieg and T illma n (reference 26) suggest 
a formula for the skin-friction coefficient, which is based on the form m 

parameter H and on the Reynolds number R^: jo 

Op » O.l23>a0“ 0 - 678H R "°* 268 

^0 0 

This equation has been applied to the data of Sehubauer and Klebanoff 
(reference 18) with the experimental values for H and RqJ The results 
are plotted in figure 13j the Ludwieg-Tillman curve falls below that of 
the hot-wire points; however, it is apparent that the two curves are 
similar, since very good correlation is obtained when the Ludwieg- 
Tillman curve is multiplied by 1.52, which is the ratio of the skin- 
friction values at the initial point (x a 17.5). 


CONCLUSIONS 

The discussion of the behavior of the turbulent boundary layer in 
adverse pressure gradient leads to the following conclusions: 

Surface roughness and stream turbulence. affect the skin friction, 
and thus should be specified in boundary- layer investigations . 

The skin friction determined according to the momentum- survey 
method cannot be considered a function of a single pressure parameter, 
as in laminar flow. A turbulent- energy production parameter is proposed 
as an additional quantity upon which the momentum skin friction may 
depend. This parameter also appears to be related to the skin friction 
as determined by the hot-wire method. 

The skin friction obtained by the momentum method for adverse pres- 
sure gradient shows an increasing downstream trend, whereas data obtained 
by the hot-wire and heat-transfer methods show a decreasing trend. 

In order to establish the correct skin-friction trend, it may be 
necessary to consider violent fluctuations and secondary flows near 
separation; it would also seem advisable to investigate thoroughly the 
accuracy of the accepted method of inferring wall friction from the 
Reynold's stress profile determined by the hot-wire method. The 
Ludwieg heat-transfer instrument shows promise as an engineering tool 
for measuring the true skin friction, and should be fully exploited. 

The momentum- survey method for the experimental determination of 
skin friction in an adverse pressure gradient is of doubtful value 
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■because this method is highly sensitive to errors in measurement and 
because secondary flaws (not contemplated in the method) probably exist 
near separation. 

Available experimental information does not permit drawing firm 
conclusions as to the true skin- friction trend in an adverse gradient. 
If, however, it is assumed! that either the momentum method or the hot- 
wire method yields the correct trend, the weight of evidence would 
appear to be in favor of the hot-wire results. 


National Advisory Committee for Aeronautics, 
Lewis Flight Propulsion Laboratory, 
Cleveland, Ohio, March 21, 1950. 
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APPEHDIX A 


SYMBOLS 


The following symbols are used in this report: 


turbulent contribution to deviation from Prandtl’s assumptions 
as given by equation (B17) 


mean flow contribution to deviation from Prandtl’s assumptions, 
as given by equation (B18) 


coefficient of shear at the wall, 



dimensionless energy transfer from mean to turbulent flow, as 
given by equation (C13) 

mean velocity-profile-shape parameter, (6*/©) 

profile-shape parameter, accounting for complete momentum 
deficiency, (5*/©-^) 

instantaneous static pressure 

Reynolds number formed with chordwise distance, (Ux/P) 

Reynolds number formed with momentum, thickness, (U©/v) 

free-stream velocity at edge of boundary layer 

first derivative of TJ along chordwise direction, (du/dx) 

second derivative of U along chordwise direction, (d^tl/dx^) 

instantaneous velocity in x- direction 

instantaneous velocity in y- direction 

instantaneous velocity in z- direction 

distance along streamwise or chordwise direction 

distance along direction perpendicular to wall 

distance along direction perpendicular to x-y plane 


2192 
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turbulent-boundary- layer thickness, extending by definition to 


6* 


e 


n 

v 

p 

T 


XX 


'xy 


yx 


point where 


u' v' 


is equal to arbitrary small value 


displacement thickness, representative of mass-flow deficiency 



momentum thickness, representative of mean- flow momentum 



deficiency. 


turbulent momentum thickness, representative of total momentum 
deficiency, as given by equation (B16) 


• (£0 


Pohlhausen parameter 

absolute- viscosity 
kinematic viscosity 
mass density 

total normal stress acting along x and on y-z plane 

total shear acting along x and on x-z plane, ^ + (-pu'v')j 

' Sv 


total shear acting along y and on y-z plane, (p — + (-pu'v 1 ) 

V dx 


'yy 


total normal stress acting along y and on x-z plane 

f 3u — 2^ 
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ijr energy thickness, representative of mean flow- energy deficiency 



^ turbulent energy thickness, representative of total energy 
deficiency, as given by equation (C12) 


Primes on small letters indicate fluctuating quantities 
Bars above symbols indicate temporal mean quantities 
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appendix b 


DERIVATION OF INTEGRAL MOMENTUM EQUATION 


The Reynolds equations of motion of an incompressible viscous 
fluid for two-dimensional steady flow with three-dimensional turbulence 
are 


du t — Su t 1 
P 


U -or — + V — + — - 5 — 

ox oy p ox 


^ a ~ • a _ 


a 

Si 


dv — Sv 1 <dp ( a^v S^v'\ B 


u — + v — + 
ox 


Sy ‘ P dy U ^x 2 + Zy 2 J~ Sx U ' V * ~ ^ V ~ Sz T W 


a 

35 


_a_ 

Sy 


a ,2 
a z 


J 


The equations of continuity are 


+ * 0 

ox oy 




> 


(B2) 


Su 1 . Sv 1 dw’ 

Sx ay + a z 

Assuming no variation of turbulent mean quantities in the Z -direction 
and assemblying the stress terms- produces 


_ du _ du 1 dp 


_ by —by 1 
u 2^ + P 


1 a 

au 

p 3x 

^ + 

1 a 

av 

p Sx 


a 

a 

— U V T 

ox oy 



1 a 

au 

(-pu 1 ) 

+ P Sy 

» ^ + (-p u v> ) 


1 a 

av # 0 

+ p¥ 

^ + (-pv^) 


(B3a) 


(B3b) 


0 
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When the following definitions are made, 


5u , , ,2\ _ - 

» s + <-p u > - T 


XX 


du 


P "Sjf + (-pu*v*) a T 


xy 


dv 

(-pu'V) a Tyx 


T 1 |^ + ("P^' 2 ) 3 T 


y y 


and the third equation of motion is neglected, the following equations 
are obtained: 



+ v 


du 1 dp 
dy + p dx 


1 S T 

p dx xx 


.id 


xy 


u 


dv 

s 


— dv 
+ v d^ + 


1 dp 
P 3jF 


id .id 
p 3 x T yx p 5y T yy 


(B4a) 

(B4h) 


When the following relations are used. 


_ du d _ _ 1 d _2 

V d^~^ U * V+ 


and 


and 


— du 

1 d -2 

Ur- a 

dx 

2^ U 

— dv 

— du 

V dy ~ 

" V dS =a 


— uu. 

' 2 t 5 T 


equations (B4) may he written: 


d — 2 , d 1 dp 1 d _ id 

SW u + ^ uv + TY: = rr T w + - -sm T 


dJ 


p dx p dx xx p T xy 


(B5a) 


1 dp 
P dy 


1 dp _ .Id 
P dx yx p dy 


T _ 

yy ds 


u v + 2 

dx 

(Bob) 
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cu 

CD 

H 

cvi 


Differentiating equation (B5b) with respect to x yields: 
1 S 2 „ . 1 ^ T yy d 2 — — . „ d v du 


. 2_ 
1 d p 


- d 2 ^ 


= + f srsy - ^s u v+ 2 ss +2T |? (B6) 


Inasmuch as the order of differentiation may "be exchanged, equation (B6) 
may he written: 


d 

3y 


^i dp^\ i s 2 T 
, P p ^x 2 ; 


S / 1 s 


- d 2 ^ - d 2 ? 


„ , W f JL IS _ 1 . — UU — U V 

r y* + ^y\p sx T yyJ + v ^2 - u Sx 2 


(B7) 


The potential flow impresses on the Boundary layer a certain pressure 

v~? 

gradient) therefore at y =» & is given hy the potential flow. The 

,ox 

is modified inside the boundary layer is 


manner in which dp 

ox 


ys8 


indicated hy equation (B7). Integrating from y = 8 to y «*• y yields 

py 


1 

p 


5 ( dp 


> 


§3^5 P " p 


dx 


1 

P 


dx 


y*5 


Inasmuch as 


Then 


dp 

dx 




i 

p 


"V 

Sy ) 


, l 
<3y = - 
p 


dp 

dx 


+ UU’ 


Therefore 


1 dp 

- UU’ + - 
P 


ny 

_d 2 _ 
v^uay- 


vj 

8 O 

6 L. 


py 


_d 2 _ 
u — g v dy 
ox 

(B8) 
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Expressing the velocity quantities as dimensionless ratios over 
the local free- stream velocity U, dividing hy U 2 , and simplifying 
yields 


~ o . _r. 

21_ + Al _d_ u ■ v dx - 2 — Txx + Txx . T xy 

17 U 2 dx n 2 + dy U U pu 2 u P U 2 dx pTJ 2 dy pU 2 


(B9a) 



py 


v 

u 



2 


ir u 
3x u 



u 


8 


py 


u ( U" v 0 U' d v ,S 2 v\ , 
U \jJ U + 2 u Sc u Sx 2 U 


(B9b) 
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dp 

'T - 

Substituting for from equation (B9b) into equation (B9a) yields: 

pU 2 


U' u 2 d u 2 


d u v U' 2 U* T xx d T xx d T xy 

^uu--u 53 "IT ^2 + dx pU 2 dy pu 2 


o IT ]yy d ]yy 

u p u 2 ' S p u 2 


P6 


y 

P5 


9 ( 22l 4 . 21 * h* + a 21 A Zz5 + 

W 2 U J pu 2 U Sx P u 2 dx 2 pu 2 


Jy 


v f U" u - U' d u d 2 u . , 

u(yr ij + 2 -u 5 5 t ^3 |air + 


dy + 


ps 


u /u" v U* d V d 2 V . 

u l u u 2 0 


(BIO) 


From the continuity requirement the following equations are obtained: 


1 21 

“ u a u 


py 


2dy + 


0 

P8 


u 


dx U 


<3y 


v 
• U 


K 21 

y=5 u 


u , , d 

u ^ 4 dS 


* dy - ® 
U ^ dx 


0 
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so that 


u 


S v • U' 


u 


Sx U + Sy U ~ U U 


Noting that 


06 


U V 

U U 


yi»5 


21 

U 


05 


u , d 
U ^ " dx 


u , , dB 

- cly + — 
U dx 


xy 


dU 2 


=3 0 


yfc*6 


and that 


06 


-2 
o u 


o ^u 2 


-o 


_d_ 

dx 


0 


u 2 , d5 

_ dy _ 

U 2 ^ 


(Bll) 


(B12) 


(B13) 


(B14) 


and integrating equation (BIO) across the Boundary layer from y =■ 0 to 
y =■ 6 yields 
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The individual terms of the complete integral momentum equation will be 

ns 


considered. The integral 


—2 
u u^ 


U 


IT 


^XX T yy \ 

„ + — 5- ^ ) dy corresponds to 


pU° pU 


0 


the momentum deficiency thickness and may be written in the following 
form: 


9 t “ 


"'S 

- 



u 

1 - H 

( x T yy 

T xx\ 

U 

1 u- 

l 1 + p * 2 

P u 2 y 

0 



j 


sy 


(bis) 


If equation (B16) is identical to the well-known laminar momentum 


thickness 9 = 


u 


T? |1 " S )dy, except for the factor (l + =£), 


1 yy 


XX ' 


u 


u 


0 


pu 2 pu 2 / 


’which represents the turbulent -momentum contribution. 

PS / 


The integral 


1 u ) 

1 ~ 1 <3y corresponds to the mass -flow 

l U / 


0 


deficiency and may be termed the displacement thickness 6* . 


The term 


'ay 

pu 2 


represents the total shear coefficient at the 


y=0 


wall, which may be composed of both laminar and turbulent contributions; 
however, at the wall itself, the shear is believed to be entirely 


laminar . This quantity 


'xy 


plT 


is termed C-p . The double integral 

y=o f o 




^5 

^y 

0 

6 


U' 

vU 2 


U' 


it \ T 


yx 


+ WW I — — + 4 : 


U’ S T yx d 2 T ; 


U 


3U 2 


U Sx 


+ 


yx 


pU Bx pU 


dy dy 


(B17) 


hereinafter termed A, represents the contribution of the anisotropic 
low-frequency turbulent field to the deviation from the Prandtl boundary- 
layer assumptions. If the turbulent field were isotropic, would 

be zero and term A also would be zero. 
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The integrals 

pSpy 


U' 

U 


Uo 




ns 


'SuV ,T S uV 

^ur 


py 


Uo, 



dy dy 


(B18) 


hereinafter termed B, represent the contribution of the mean flow field 
to the deviation from the Prandtl boundary- layer assumptions. 


In conclusion, the complete integral momentum equation may be 
witten as follows: 


d9 t 

dx 


+ ~ (20. +5*) « C, 


U 


+ A + B 


Por comparison the Karmen momentum equation is 


d9 U* 
dx + U 


( 20 + 6 *) 
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APPENDIX C 

DERIVATION OP INTEGRAL ENERGY EQUATION 

If the Navier -Stokes equations are multiplied by the correspond- 
ing instantaneous velocities and time averages are taken after the 
instantaneous-energy quantities are formed, equations describing the 
average energy balance of a two-dimensional mean flow with three- 
dimensional turbulence are obtained: 


n 2 Bu chi , chi 

U— tU^ + uVv + WUr- 
ot ox oy dz 


p OX 


/c)^u B^ B'nA 

W By 2 Bz 2 / I 


vll + uvl^+v^il+vw 


Bv 

Bx 


¥ 


Bv 

33 


z ^ + w (¥i + ^!z + bM 

P °Y \Bx 2 By 2 Bz 2 / 


> (Cl) 


Bw 


w^ + wu^ + wv^ + w 2 |Z=-Z^ + w + f^W + 


\Bx 2 By 2 Bz 2 / -S 


Equations (Cl) can be rewritten as follows: 


1 Bu 2 u Bu 2 + v Bu 2 w Bu 2 

2 Bt 2 Bx 2 By 2 Bz 


1 Bv 2 u Bv 2 v Bv 2 w Bv 2 

2 Bt 2 Bx 2 By 2 Bz 


E |E + uu 

p Bx 

v Bp , 2 

— + uvyv 

p oy 


1 Bw 2 u Bw 2 v Bw 2 w Bw 2 _ __ _ _ 

2 Bt 2 Bx 2 By 2 3z p 3z 




> (C2) 


= - Z^E+ywy 2 ^ j 


Adding the three equations yields : 


1 (u2+v 2 +w 2 ) + | ^ (uW+w 2 ) + f (u 2 +v 2 +w 2 ) + 


2 31 


2 3^ 


1 JL (u 2 +v2+w 2 ) = 

2 Bz 


“pfli +v Bf + w l^) + u ( u 7 2u + v Vv + ^ v 2 w) 

(C3) 
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and 


The following relations may now be written: 

u = u + u' 
v = v + v* 
w = w + w' 

P = P + P' 

u 2 + v 2 + w 2 = q 2 
u' 2 + v >2 jj. w i2 _ g»2 


It is then assumed that w = 0 for two-dimensional mean flow; thus 

u 2 + v 2 = q 2 

u' 2 + v' 2 + w' 2 = q' 2 

and 

q 2 = (u+u 1 ) 2 + (v+v') 2 + w' 2 


q 2 = q 2 + 2 (tru'+w') + q 


,2 


The left-hand side of the equation (C3) thus becomes: 


1 

2 eft 


£q 2 + 2(uu'+w') + q' 2 J + t _ A . A. |q 2 + afuu'+w') + q ,2 J + 

— — — A j^j 2 4. 2(uu'+w') + q'^ + ^ Jjl 2 + 2 (uu'+tv i ) + q ,2 j 

Considering the time derivative, analyzing term by term, and taking 
time averages yield: 


1 dq 2 _ 


2 3t . 


= 0 
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because the mean motion is steady, and 


because the time average of uu’ is zero. 
Similarly, 



and 


1 Sq' 2 

2 St 


1 IP> - 


0 


because, by definition of the mean, the averaged turbulent quantity 
is assumed to be independent of time. Analysis of the terms contain- 
ing x-derivatives yields 


u dl 2 _ u Sq 2 
2 Sx 2 Sx 


because u' = 0, 


— Su u’ _ S /="Tn „ 

u — ^ = u tMuu') = 0 

ox ox 


_ Sv v' _ S . 

u — - — = u v-^w’) = 0 

Sx OX 


because v* = 0, and 


u Sq' 2 u S . 

zsr = 5 ^ > 
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The quantity pu ^ 




is the net gain of kinetic- energy 


transport in the flow direction. 
Now: 


u * do 2 u * So 2 

— = o because H* = 0 by definition 

2 dx 2 dx 


. , d(uu') , / , du , du'\ — ; — 7 du , u d(u'u') 

dx \ OX dx / bx 2 bx 


u 


, d(w’) t f , bv , — dv’^ — ; — 7 

= u ( v s + T srj= u ^ 


1 — r dv _ fb(u'v') , du 1 ^ 

v si + v \“ L s~ " v srv 


i 1 2 

The quantity - is a turbtilent energy convection term and 

is thus obtained for the x-derivative terms: 

— d / q^ + q'^ , — r — r du — — : bv u dfu’u 1 ) 

u i ^ — l+u'u 1 + u'v' Xi + h. + 

ox \ 2 / bx Si 2 bx 


d(u' v 1 ) _ , du 1 


^i 


dx 


u ' dq ' 2 

T-sr 


Analysis of the y-derivative terms yields: 


▼ dq 2 _ v dq 2 
2 dy 2 dy 


v = v (uu ’ ) = 0 

dy dy 


because H' = 0 by definition, and 


dw' — d /— , v 

V -5T = T 3? (TT 5 = 0 


because 


v* = 0 by definition 
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v q' 2 _ v Bq' 2 
2 y 2 By 


— B ^q 2 + q' 2 ^ 

The term pv ^ y ^ y is the net gain of kinetic energy of the 

mean flow in the y direction: 


v 1 Bq^ _ V Bq 2 _ q 
2 By 2 By 


because v 1 = 0 by definition 


, B(uu') — : — : Bu , Bu 1 — : — ■, du — 

v' — - = v'u' — + v'u - — = u'v' — + u 

By By By y 


B(u'v') _ Bv' 

~ 57 ~ W 


, B(w') — — ; Bv , v B(v’v') 

T - T ' T ' 3F + 2 


v 1 Bq* 2 _ Bq' 2 
2 By 2 By 


v 1 Bq 1 2 

The term — — — is a turbulent energy convection term. Thus, for 
2 By 

the y derivative terms, there is obtained: 


v -r- —■ — + u'v' + v-v- + u 


- B 
‘5y 


4- iT 

By 


d(u'v , _) _ u , V 


dy 


v B Cv’v') v]_ Bq' 2 

2 By 2 By 


Analysis of the terms containing the z derivatives yields: 


o=0 

2 Bz 2 


because v' = 0 


and = 0, 

oz 
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■w* 


B(uu') — — ; 

— - = w'u* 
Bz 


Bu , Bu 
Bz U Bz 


u 


B^u'w* ) 
Bz 



because 





M 


w 




— = w'v' 


Bv 

c)z 


+ w’v 


Bv' 
3^ 



because 


— = 0, and 
Bz 


v 1 Bq,' 2 _ Bq' 2 
2 Bz 2 Bz 


Thus, for the z-derivative terms, 


r 1 


r 1 

BCu’w*) 

- u 1 

+ V 

B(v , -W*) 

Bw' 

Bz 

Bz 

Bz 

" v ST 


v' Bq' 2 
+ T "5T 


When the pressure terms in equation (C3) are considered and moved to 
the left side of the equation. 


1 

P 


<^ ,) (ii:& : ) +(?w,) fe + £ L ) +w, £ 


*v=r 

because = 0, and 


i( S ^ + u|El + T^ + 7^- + u’|E +T ’p + u'3P: T . + 

p \ ox Bx By By Bx By 3y 


w 


, Bp 

Bz 


-) 


Taking time averages and noting that 5' = 0 yield . 

l /_ Bp" _ Bp , Bp’ , Bp 1 , Bp’ ) 
p \ dx By Bx By Bz ' 
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Thus, the complete left side of the equation, including the pressure 
terms, can be written as follows: 



(a) 


a 


o 


u 1 3q |2 Bq' 2 w^_ 3q. t2 u_J_ Bp 1 y]_ Bp* w]_ Bp 1 /, v 

2 5x 2 3y 2 Sz p Bx p By p 5z 


u'u' + u'y' ^ + u'y' + v'v' + 
ox ox By By 


(c) 


r . . 1 


r 

1 


u B(u'u') . Y B(v'y') . — 13 

(u'v) T 

Bu 1 

y t i 

BCu'v*) Bv 1 


2 Bx 2 

dy L 

Bx 

Bx 

*r u 

_ dy 

dy. 

+ 

u 

d(u'w’) , 

dw ' +v 

B(v'w 

li - 

. Bw’ 



. az 

Bz J 

Bz 



j 

(a) 


The next problem is to reduce the expressions to a more compact 
and general form: 

The group of terms listed as (d) can be rewritten as follows: 


1 

11 2 


BCu'u*) dOu'v* ) BQu'w 1 ) _ , /By 1 3w'\ 

Bx By Bz \3y Bz / 


1 5(-y t y l ) B^'y') 5(w , w t ) , Au' Bw 1 \ 

2 By + Bx Bz V \3x + Bz / 


u 


dCu'u') , ^'v') , ^fu'w*) _ u , SuJ. _ i,' ^Xi _ „i Sw* 

Bx By Bz Bx By Bz 


BCy'y’) B(u’v') BCy’v') . Bv 1 . Bu' . Bw' 

— ^ + v <■ + \ ' - Y* -s - v' -s V 1 3 

By 6x Bz By Bx Bz 
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because 


1 d(u'u') _ u , + u t d^ _ u , d^J_ _ SCu'u') _ u , ctoi 


35" ^ u 35" “ u 35 3F 


dx 


and 


1 = v* |ll + V |ll - v . _ v , 3v 


Sy 3T " V 3T = ~^y" 


ST 


3u 1 3 v * 3w 1 

Inasmuch, as, by continuity, -^— + = 0, the following 

equation is obtained: 


•a 


d(u'u’) . dfu'v') . d(u'w*) ,/du' . 3v' dw' ^ 

ar^ + + - u \sr + sr + srV] 

/du' dv 1 dw 1 ^ 

\dx dy dz /_ 


3(u t v' ) d(v' v 1 ) d(v'w' ) 

dx dy dz 


u 


d(u'u') d(u'v') d(u'w’) 

4- Tr 

d(u'v') , d(v’v') , d(v'w') 

dx + dy + dz 

T V 

[& + a y + a z J 


(*) 


If the group of terms in (c) is conib ined with the new group of terms 
in (d), then 


uhT 7 + u ^L+v 5 ^ ,v> ) + 

ox ox ox ox 


— r du — d(u'v') —. — r dv _ d(v'v') , 

oy dy dy dy 


■ p —i du — d(u’w’) ■ 
dz dz 


1 — r du 

I 1 1.T • 


C*U 


since u'w' ^ = 0 because g— = 0 


— ; — , dv , _ d(v’w’) 
dz dz 


dv 


du 


since v’w' = 0 because = 0 
dz dz ■ 
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or adding 

4L(u u 'u * + t 11 1 v' ) + ^_(u u'v' + v "v 1 v 1 ) + ^_(u u'v' + v v'-w') 
dx dy dz 

The group of terms listed as (b) may be -written as follows: 


u' dq'^ u' dp' v* dq'^ v* dp 1 w' dq 1 ^ w' dp ' 

*2 dx + p dx + 2 dy + P dy dz p dz 



The group of terms listed as (a) can be rewritten to yield: 
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By continuity, 


a o; thus terms (a) reduce to: 

dx dy 



It is now possible to give the complete left side of the energy 
equation in its final form: 



The terms on the right side of the equation are now considered: 

v y 2 u + v y 2 v + w l^w^ U y 2 ^ 

■— — n-(u+U ’ ) + — — U+U ' ) + — — ^-(u+U 1 ) 

^ a? J 

r-, tS f d^u d^u d^u' d 2 ^ d'hi'^ 


= (u+u’) 



dV 

\ _/a 2 u* , dV 

, d 2 ^ ’ 

\ , /d^u 

d^N 

(i? + 

dy 2 ' 

) +u L 2 + , 2 

Vox dy 

+ 2 
dz / 

) +U (d? 

dy V 


+ 


u 1 (&y± + 

\dx 2 dy 2 dz 2 / 
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When time averages are taken, the following terms are developed: 



Similarly, the other terms can be developed: 


,/d 2 v , b 2 v\ , ,/B 2 v* , d 2 v'\ 5 2 v'^ 
W 2 by 2 ) W 2 Sy 2 dz 2 ) 


and 




dV , d 2 w'\ 
by 2 bz 2 / 


Thus, the expression on the right side will be 


v ^u y 2 !! + v V 2 ^ + u’ + v* V 2 v* + w' V 2 w' ) 


(C5) 


This term represents the dissipation of flow energy into heat. [The 
dissipation mechanism is laminar, as can be seen from the form of the 
expression. 


Writing the three basic 
three-dimensional turbulence 


equations for a two-dimensional flow with 
yields 


Continuity equations: 


cfu bv 


0 


du 1 3v' bw' 

Sx by Sz 



(C6) 



NACA TN 2431 


39 


Momentum and continuity equations: 
cftr _ Su 1 Sp , S^u\ 5/— — 

E s-5 + fs = 'b t ^'S (uu) 


'dx by 

_ Sv _ 5 V 1 Sp (cPv d /— ; — ;■« 

u s +T s + Pbj=' , (^ + i 7 ;- 3 ^“ , ' r ’) 




k (r ' V,) ' ^ (V ' W,) 


> 


^(u'V) + + |^v'v') = 0 


az 


J 

(C7) 


Energy and continuity equation: 
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_3_ 

by 


u (^-— | - h + U< ^2~ + p“) + uu'u' + V u'v’ 
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l^2~ + ^p”) + u u’w’ + v u'v' 

= u ^ y^u + v v^v + u* v^u* + v* V 2 v l + w' V 2 w') 


(C8) 


Because momentum) energy, and continuity conditions must be 
simultaneously satisfied, the momentum conditions can be inserted in 
the energy equation and equations of both momentum and energy, (and 
also continuity) can be obtained. 


In the energy equation the substitutions 


“(§ + 0)-s (u ' u ' 


1 dp _ _ ;r Bu y du 

p 5E " 3x 

and 




3^ 


|.(T'V') - |(v¥) 
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directly yield 



Either the energy equation and the momentum equations or the tiro 
energy and momentum equations can therefore he used. 


Continuity, momentum, and energy equations: 
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and 



For the purposes of this report , equation (CIO) is integrated 
across the boundary-layer section from y = 0 to y. =» 8. This 
equation of momentum and energy has been chosen because it deals 
primarily with mean quantities. The viscous dissipation on the right 
side of equation (CIO) can be neglected in comparison with the 
turbulent- energy production. On the assumptions that 3p/Sy is 
actually negligible, that 


uIe 2 HI v ' 2 d / u \?' 2 , ? 2 a M 

u 3 p dx " u U " U iJ?" + S\jjy 

and that u 2 +v 2 = u 2 , equation (CIO) is integrated to the following 
form: 



The right side of the preceding equation is twice the integrated 
turbulent- energy-production terms from equation (C9); the factor 2 is 
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introduced "because on the left side there is twice the kinetic- energy- 
loss . 


’'5 
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1 - — f 

0 u'2-v ,2> \ 
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u 2 \ 

u2 /J 


dy = i|r t 


(C12) 


and 


r-PS 
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uO 



dy + 


5 u 3 /u ,2 -v' 2 a 
5x u 3 \ u 2 / 7 


u 3 /u-MV _ „ 

Ij3 \u 2 )^ 


(C13) 


give the integral momentum and energy equation as 

gt + 3 == - 2E (CIA) 

or 

a§- *t u 3 = - zuSe 

Physically, the preceding equation should he written 

s (§ pDS *t) - - pu3e 



13 

The loss of kinetic energy g- pU ^ is equal to the total turbulent- 
energy production over the houndary-layer volume plT^E dx. 
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Figure 2. - Mean velocity profiles at same station for various free-atream tur- 
bulence levels. Data from reference 20. 
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(b) Data from reference 18. 

Skin -friction parameter under adverse pressure gradient. (Skin 
friction computed by momentum method.) 


Figure 3. - Concluded. 



Skin-friotlon coeffioient, Cf 
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Figure 4. - Comparative momentum and hot-wire skln-frioti on coefficients under adverse 
pressure gradient. Data from reference 18. 
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Figure 7. - Relation of skin friction and pressure parameter for two different 
free- stream turbulence levels under adverse pressure gradients. Data frcra 
reference 21. (Skin friction computed by momentum method.) 
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Turbulent-energy parameter, E 

Figure 9. - Relation of skin-friction parameter and turbulent -energy parameter 
for constant -values of pressure parameter. Data from references 18 and 25. 
(Skin friction computed by momentum method.) 
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0 ,-t .8 1.2 1.6 2.0 2 . 4 , 2.8 3.2 3.6 

Distance from wall, y, in. 


(a) Chordwise distance x, 17.5 feet. 

Figure 11. - Mean and turbulent velocity profiles for computation of the laminar and turbulent 
momentum thicknesses based on data from reference 18. 



Turbulent velocity ratio 
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(b) Chordwlse distance x, 20.0 feet. 

Figure 11. - Continued. Mean and turbulent velocity profiles for computation of the la m i n a r and 
turbulent momentum thicknesses based on data from reference 18. 
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(o) Chordwise distance x, 21.0 feet. 

Figure 11* - Continued. Mean and turbulent velocity profiles for computation of the laminar and 
turbulent momentum thicknesses based on data from reference 18. 
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(a) Without turbulence grid. 

Figure 12. - Mean and turbulent velocity profiles for computation of the lamin ar and 
turbulent momentum thicknesses based on data from reference 20. 
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Distance from wall, y, mm 
(b) With turbulence grid. 


Figure 12. - Concluded. Mean and turbulent velocity profiles for computation of the 
laminar and turbulent momentum thicknesses based on data from reference 20. 








